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V SUMMARY 

A finite  element  method  for  analysing  the  behaviour  of  cracks  in  metal  sheets,  which 
are  patched  with  an  overlay  of  composite  material,  is  described.  The  analysis  includes  the 
separate  responses  of  the  sheet,  composite  patch  and  the  adhesive.  Debonding  of  the 
adhesive  is  permitted  and  the  stress  intensity  factor  at  the  crack  tip  is  obtained  as  well  as 
the  stress  distribution  in  the  patch. 
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ABSTRACT 

A finite  element  method  for  analysing  the  behaviour  of  cracks  in  metal  sheets,  which 
are  patched  with  an  overlay  of  composite  material,  is  described.  The  analysis  includes  the 
separate  responses  of  the  sheet,  composite  patch  and  the  adhesive.  Debonding  of  the 
adhesive  is  permitted  and  the  stress  intensity  factor  at  the  crack  tip  is  obtained  as  well  as 
the  stress  distribution  in  the  patch. 
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Direct  stresses  in  the  patch. 

Transverse  shear  stresses. 

Elemental  displacement  vector. 

Angle  measured  from  local  axis  system  to  the  material  symmetry  axis. 

Distribution  function  of  the  transverse  shear  stresses  within  an  element  of  the 
adhesive. 

Transverse  shear  strains. 

Strains  in  the  patch. 

In-plane  coordinates  and  the  coordinate  in  the  thickness  direction  respectively. 
Thicknesses  of  the  sheet,  adhesive,  and  patch  respectively. 

Poisson’s  ratio  for  the  sheet. 

Shear  moduli  of  the  sheet  and  adhesive  respectively. 

Transverse  shear  moduli  of  the  patch  in  the  xz  and  yz  planes  respectively. 
Displacements  at  the  mid  surface  of  the  patch. 

Displacements  at  the  mid  surface  of  the  sheet. 

Strain  energy  of  the  adhesive. 

Mode  1 stress  intensity  factor. 

Element  stiffness  matrix  for  the  adhesive. 

Element  stiffness  matrix  for  the  patch. 

Element  stiffness  matrix  for  the  adhesive-patch  pair. 
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1.  INTRODUCTION 

A bonded  overlay  of  high  strength  composite  material  offers  an  efficient  method  for  reducing 
the  stress  intensity  factor,  and  hence  the  stress  field  at  the  tip  of  a crack.  This  reduction  in  the 
stress  intensity  factor  subsequently  retards  crack  growth.  The  primary  objective  of  the  work 
reported  here  is  to  develop  a method  of  analysis  for  a crack  which  has  been  patched.  The  develop- 
ment is  based  upon  the  finite  element  method  for  analysing  cracked  structures  as  previously 
described1,2  by  the  authors.  The  present  analysis  is  an  extension  of  this  previous  work  to  patched 
cracks,  and  it  allows  for  separate  responses  of  the  adhesive  and  the  composite  patch  as  well  as 
the  debonding  of  the  patch  from  the  sheet. 

This  report  forms  the  first  section  of  a detailed  investigation  into  crack  patching  and,  as 
indicated  above,  its  main  purpose  is  to  develop  the  theoretical  tools  for  the  investigation.  A 
detailed  study  into  the  effects  of  changing  the  geometry  of  the  patch  as  well  as  changing  its 
material  properties  is  to  be  reported  in  a subsequent  paper. 


2.  METHOD  OF  ANALYSIS 


Let  us  begin  by  considering  a thin  composite  patch  which  is  bonded  to  a thin  clastic  sheet. 
The  x and  y axes  are  taken  in  a plane  parallel  to  the  mid-surface  of  the  sheet  while  the  z-axis  is 
in  the  thickness  direction.  The  sheet  is  subsequently  acted  upon  by  a system  of  external  loads 
acting  at  its  edges  and  with  lines  of  action  lying  in  its  plane.  When  the  patch  is  not  present  this 
problem  is  usually  tackled  under  the  assumptions  of  generalized  plane  stress.  However,  when  the 
patch  is  present  shear  stresses  will  be  developed  in  the  adhesive  bond  and  it  is  reasonable  to 
assume  that  these  will  be  continuous  across  the  adhesive-sheet  interface  as  well  as  across  the 
adhesive-patch  interface.  Furthermore  these  shear  stresses  tzz,  and  tvz  are  zero  at  a free  surface 
or  at  a plane  of  symmetry,  and  it  is  reasonable  to  assume  that,  since  the  patch  and  the  sheet  are 
thin,  these  stresses  vary  linearly  with  thickness  in  the  patch  and  the  sheet.  With  these  assumptions 
the  distribution  of  the  shear  stresses  rzz  and  t„z  is  found  to  be: 


tX2,  Tyj  — Tnz,  Tgy  for  Z ^ 
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where  t,z , r,y  are  the  shear  stresses  in  the  adhesive  and  are  assumed  constant  through  the 
thickness  of  the  adhesive.  The  2 — 0 plane  is  taken  at  the  mid-surface  of  the  sheet  for  a doubly 
reinforced  sheet  and  at  the  lower  surface  of  the  sheet  for  a singly  reinforced  sheet.  Here  l0,  la 
are  the  thicknesses  of  the  patch  and  adhesive  respectively  while  t,  is  the  thickness  of  the  sheet  in 
the  doubly  reinforced  case  and  is  twice  the  sheet  thickness  in  the  singly  reinforced  case.  This 
stress  distribution  is  shown  in  Figure  I. 

In  what  follows,  we  will  confine  our  attention  to  the  symmetrically  patched  sheet  but  we 
note  that  the  results  can  also  be  applied,  with  only  minor  modifications,  to  a singly  reinforced 
sheet  of  thickness  r„/2. 
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FIG.  1 DISTRIBUTION  OF  SHEAR  STRESS  THROUGH 
THICKNESS  OF  PATCHED  SHEET 


Let  us  now  consider  the  x-y  axes  to  be  at  some  angle  8 to  the  axes  of  orthotropy,  denoted 
by  x for  the  composite  patch  (see  figure  2). 


FIG.  2 AXIS  SYSTEM  IN  PATCH 

In  the  patch  material  the  transverse  shear  stresses  tI7,  tv.z  are  related  to  the  shear  strains 
Yx-t , yyi  through  the  formulae 

tx-t  = C'l3  Yz’t  (2.1) 

*yx  = G' 23  Yvi  (2.2) 

where  G'u  and  G' 23  are  the  interlaminar  shear  moduli  for  the  patch,  while  the  shear  stresses 
rxt,  Twt  are  related  to  rz.f,  ry.t  as  follows: 

t xt  = t x’z  cos  8 — Tyi;  sin  8 

Ty,  = t x't  sin  8 + Ty.t  cos  8 


2 


(2.3) 

(2.4) 


Hence  substituting  equations  (2.1)  and  (2.2)  into  (2.3)  and  (2.4)  we  obtain 
r xi  — G'i 3 yz-i  cos  # — G' 23  Yv2  sin  0 
tvi  = G'13  yx-z  sin  8 + G'23  y% rz  cos  8 

However, 

yx’z  = yxz  cos  8 4-  yyz  sin  8 
yt rz  = —Yzz  sin  8 + yvz  cos  8 

so  that  on  substituting  for  yx-z  and  yv-z  into  (2.5)  and  (2.6)  we  finally  obtain 

txz  = G'  13  ( yxz  cos "6  + yyz  sin  8 cos  8)  — G'23  (yyz  cos  8 sin  8 yzz  sin2#) 
Tyz  = G'13  (yxz  cos  8 sin  0 4 yyz  sin2#)  4-  G'23  {yyz  cos2#  yxz  sin  8 cos  8) 

Now, 


Yxz  = 


iu 
iz  ’ 


iV 


Yvz  = . 


(2.5) 

(2.6) 

(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 


where  u and  v are  the  horizontal  components  of  displacement  in  the  .v  and  y directions  respectively. 

Substituting  for  txz,  and  tvz , as  given  by  equation  (2.1),  and  for  yzz  and  yvz,  as  given  by 
equation  (2.1 1),  yields 


where 


_ iu  iV 

Az)r.„  = -fm  + M8) 

is  is 

hU  iv 

f(z)  T.X  = - M8)  + /.(#) 

vZ  vZ 

/i(#)  = (G'13  — G'23)  sin  8 cos  8 
f-z(8)  = G'23  cos2#  4-  G'13  sin2# 
fz(8)  — (G'13  cos2#  (-  G' 23  sin2#) 


/(’) 


- + U f t0  - z) 


(2.12) 

(2.13) 

(2.14) 


and  where,  since  we  are  examining  the  patch  behaviour  ta  4-  h/2  ^ z ^ ta  4-  t0  4-  /*/2.  If  we 
multiply  equation  (2.12)  by  /!(#)  and  equation  (2.13)  by  fi(8)  and  subtract  we  obtain 


(fm  -/2(#)/3(#)]  =/(z)  [/,(#)  T,v  -f2(8)  rtz] 

iz 


(2.15) 


Similarly  upon  multiplying  equation  (2.12)  by  f*(8)  and  (2.13)  by  fi(8)  and  subtracting  we  obtain 


[/l2  -f2f3)lV  =/(*)  [fszx  -f3rsV] 
OZ 


(2.16) 


Recalling  that  /(z)  = (t„/2  + ta  + to  — z)/to  and  integrating  equation  (2.15)  with  respect  to  z 
yields 


. . Z\(/lT«0  —fiTtx)  , . . 


(2.17) 


where  h(x,y)  is  an  arbitrary  function  of  x and  y.  If  we  now  denote  the  movements  of  the  mid- 
surface of  the  patch  by  ug,  vB,  i.e.. 


u0  = u(x , y,  z)  at  z = 4-  ta  + 


\ f f « 

t’o  = t<jr,  y,  z)  at  z = r 4 4-  t 


(2.18) 
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we  then  see  that 


.«(£+/.  + 2)  (4+2 


a . 3fo\  (f\r,y  — fiTgz)  . 

2 4 / C/i*-/*/8)la  (*’ >) 


so  that 


and  hence 


/|(X’,  >>)  = Wo 


(/iT«v  — hr, z) 
(/i2  — hh)io 


. (hT’u  — hT>t)  [ / /v  , , 

“ " f tff  -MK  Kj  1 '*  1 ' 


K?  +'•  + '•-  0 - (j  + '•  + '2')('4'  > j + ’4-)}  (2  21) 


Similarly  we  find  that 

, f/i-'M  - h 


Hi Ki ' u 1,0  - 2)  ~ (2 1 ,a + 2X4 1 2 1 4°)}  (2-22) 


“ (/.2-/2/3)/o  l \2  ° V \2  - 2/\4  ’ 2 1 4 /) 

The  w and  v displacements  at  the  patch-adhesive  interface  c i,/2  -f  ta  are  therefore  found  to  be 
. (hTsy  hrsx)  f ( t,  \(  t,  ta  \ ( t,  1 0\/ 1,  ta  3/o\l  «... 


u = u0  \ 


V = v„  -t- 


ts  f o\(  fs  fa  3/o\|  _ _ 

2 +/a+  2/(4  + 2+j)  ( ) 


(f\r,y  hT*z)  t>  . , \(  **  ta  \ / t$  I o\f  t,  ta 

w-mio  1(2  f 'aJ{4  f 2 + ( 2 ’ 2/(4 + 2 

(f\r,x  hT»v)  \(  t,  \( I’  *a  \ 1 1,  to\(  t,  ta 

</,*-/>/>)'.  lU  ‘A*  1 ) v 2 1 '■  : 2JU  1 2 


(f\r,x  —j 


+ t)}  (2-24) 


Let  us  now  turn  our  attention  to  the  shear  ‘"-ess,  and  the  displacements,  in  the  sheet  which 
we  consider  to  be  isotropic.  In  this  case 


iz  G,  T“ 


iz  G.  ) 

where  G,  is  the  shear  modulus  of  the  sheet.  Substituting  for  the  shear  stresses  rxz<  Tyz  as  given 
by  equations  (2.1)  and  integrating  with  respect  to  z we  obtain 

.2 

“=  ” r,x  + g(x,y)  (2.26) 

G,t, 

where  g(x,  y)  is  an  arbitrary  function  of  x,  y.  Denoting  the  displacements  of  the  mid-surface  of 
the  sheet  (i.e.  z — 0)  by  u„,  v,  we  obtain 


u,  = g(x,  y) 


so  that 


and  similarly 


w = u,  + 


V — V,  + 


The  maximum  value  of  the  displacements  occur  at  the  sheet  adhesive  interface,  z — t,j 2,  and 
have  the  values 

1 ,s  T*x  im 

w = w,  + ---  (2.30) 

4G, 


V 


(2.31) 


v = 


i’«  + 


h T,V 

4 G, 


We  now  k ow  the  displacements  at  the  interfaces  of  the  patch-adhesive  and  sheet-adhesive, 
and  are  in  a position  to  evaluate  the  shear  strain  and  shear  stress  in  the  adhesive.  Let  us  begin 
by  assuming  an  isotropic  adhesive  in  which,  as  mentioned  earlier,  the  shear  stresses  tzz  and  tvz 
take  a constant  value  across  the  thickness  of  the  adhesive.  Since  the  adhesive  is  isotropic  we  know 
that 


and 


Tg*  lu  u(x,  y,  1,12  4-  tg)  — u(x,  y,  t,/ 2) 
G a Sz  la 


(2.32) 


Tty  __  iv  _ v(x,  y,  l»l 2 + tg)  — v(x,  y,  I,/ 2) 
Ga  t)Z  la 


(2.33) 


where  Ga  is  the  shear  modulus  of  the  adhesive,  and  where  u(x,  y,  I s/2),  v(x,  y,  t,/2)  are  the  values 
of  the  displacement  at  the  sheet-adhesive  interface  as  given  by  equations  (2.30)  and  (2.31)  while 
u(x,  y,  tsl 2 + ra),  v(x,  y,  t ,/2  + ta)  are  the  displacements  at  the  adhesive-patch  interface,  as 
given  by  equations  (2.23)  and  (2.24).  Substituting  these  values  into  equations  (2.32)  and  (2.33) 
and  collecting  terms  we  finally  obtain 


and 


JY  ta  . It  \ Stgft  1 . ^tofiTgy 

u„-ut  = rm  + 4GJ  - 8(/l2  fih)\  + 8(/l2  _/2/,) 


V0  — Vs  — r sy 


/«  \ 3/0/3  1 

4 G,)  ~ 8 (/,*  AA)j 


2>tofirsz 

8(/i2 -/*/*) 


These  can  be  solved  for  and  t„v,  giving 


(2.34) 


(2.35) 


T«I  = 


Tsy  — 


fa(u0  ~ Us) 

AAo  — Vs) 

(2.36) 

eg 
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(AA  -A2) 

MU0  - Us) 

AO’o  — n„) 

(2.37) 
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>> 

tc 

(AA  A2) 

where  for  convenience  we  have  denoted 


3/a/i 

W*  ~A2) 

! a Is  3/„/l, 

Ga  + 40,  ~ 8(/,2  -/2/3) 

/a  t s Hof, 

Ga  ; 40’,  ~ 8(/,2  -/2/3) 


(2.38) 


(Note  that  for  the  case  of  a singly  reinforced  sheet  the  term  /„/4G,  in  the  expressions  for  A and  /« 
have  to  be  replaced  by  3/s/8G,). 

We  thus  see  that  for  a composite  patch  bonded  to  an  elastic  isotropic  sheet  the  shear  stress 
developed  in  the  adhesive  depends  upon  both  of  the  relative  displacements  (u0  — u»)  and  (r0  — u») 
between  the  patch  and  the  sheet.  In  the  case  when  G'i3  = G'23  (=  Go),  which  only  occurs  when 
the  patch  is  transversely  isotropic,  then  the  analysis  simplifies  considerably.  In  this  case 


fiW  = 0 

MO)  = Go 
fa(0)  = Go 

/4  = 0 
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(2.39) 


so  that 


/»=/• 


t a 1 1 3/o 
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(2.40) 


(2.41) 


In  this  case  the  shear  stress  t„x  depends  only  on  the  quantity  u0  — u,,  while  r,„  depends  only  on 
t'o  — v,.  Such  a behaviour  is  known  to  occur  when  the  patch  is  isotropic.  However  in  analysing 
the  effect  of  stiffeners  on  crack  growth  several  investigators — e.g.  Arin3-4— have  assumed  the 
relationship  to  have  the  form 

t*z  ~ ■ ' (wo  — w*)  (2.42) 


(2.43) 


In  so  doing  they  have  effectively  neglected  the  contributions  of  the  terms  r*/4G„  + 3/0/8G0.  An  1 

estimate  of  the  error  involved  in  their  approximation  may  be  obtained  by  considering  a Boron- 
epoxy  laminate  of  thickness  0-508  mm  bonded  to  a 2-286  mm  thick  aluminium  sheet  the 
thickness  of  the  adhesive  being  0-1016  mm.  The  material  parameters  for  the  Boron-epoxy  are 

Gi3  = 7-24  X 103  MPa  while  the  shear  modulus  of  the  adhesive  and  the  aluminium  are  J 

9-65  X 102  MPa  and  27-3  X 103  MPa  respectively.  These  values  correspond  to  a fairly  typical 
patch  configuration.  This  gives  the  value  of  G0/ta  as  9-50  X I03  MPa/mm  while  the  ratio 


takes  on  the  value  6-56  x 103  MPa/mm,  where  for  the  sake  of  comparison 


we  have  approximated  G 23  by  G13  thus  assuming  the  patch  to  be  transversely  isotropic.  The  error 
introduced  by  using  the  simple  formulae  (2.42)  and  (2.43)  for  the  shear  stresses  is  therefore  45",,. 
This  casts  doubts  on  the  results  and  conclusions  of  those  investigations,  such  as  those  of  Arin3-4, 
which  utilize  these  simple  formulae.  Indeed  any  analysis  based  upon  equations  (2.42)  and  (2.43) 
will  significantly  overestimate  the  load  transmitted  to  the  patch  and  consequently  underestimate 
the  stress  intensity  factor. 

It  is  interesting  to  note  that  a relationship  similar  to  that  given  by  equations  (2.40)  and  (2.41 ) 
is  also  presented  by  Mitchell  et  al.5  although  no  exact  derivation  is  given  there. 

It  is  important  to  note  that,  unlike  references  [3,  4.  5],  (he  present  analysis  clearly  shows 
that  the  shear  stress  in  the  adhesive  t,z  (and  similarly  rKy ) is  influenced  by  the  relative  motions 
of  the  sheet  and  the  patch  in  both  the  x and  y directions. 

Let  us  now  turn  our  attention  to  the  case  of  a patch  which  has  G13  a G23  ( = Go)  and  for 
which  the  shear  stresses  in  the  adhesive  are  given  by  equations  (2.40)  and  (2.41).  We  may  now 
substitute  these  expressions  for  the  shear  stresses  into  equations  (2.30)  and  (2.31)  and  obtain  the 
* following  expressions  for  displacements  in  the  sheet,  viz. 


3.  K.  Arin 

4.  K.  Arin 


5.  R.  A.  Mitchell, 

R.  M.  Wooley  and 
D.  J.  Chivirut 


A plate  with  a crack  stiffened  by  a partially  debonded  Stringer.  Eng. 
Frac.  Mech.,  vol.  6,  pp.  133-140,  1974. 

A note  on  the  effect  of  lateral  bending  stiffness  of  stringers  attached 
to  a plate  with  a crack.  Eng.  Fract.  Mech.  vol.  7,  no.  1,  pp.  173  179, 
1975. 

Analysis  of  composite-reinforced  cutouts  and  cracks.  A.I.A.A.,  13,  6, 
744-749  (1975). 
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u(x,  y,  z)  = ua  + z'lK,(u0  — Ug)/t,G „ 
v(x,  y,  z)  — v,  -f  z2Ks(v0  — v,)/ttGs 
where  for  convenience  we  have  written 

3t0 


K,  = 1 

The  mean  stress  in  the  sheet  is  thus  given  by 
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and  similarly 
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(2.44) 

(2.45) 

(2.46) 


(2.47) 

(2.48) 

(2.49) 


where  v,  and  Eg  are  the  Poisson’s  ratio  and  the  Young’s  modulus  for  the  sheet,  t8Z!,,  o-si,  cr„v 
are  the  stresses  in  the  sheet  as  predicted  by  the  theory  of  plane  stress  (i.e.  when  the  transverse 
shear  stresses  are  neglected)  and  y0xy,  t ox , toy  are  the  strains  in  the  patch. 

In  the  case  of  a cracked  sheet  the  strains  in  the  patch  are  finite  while  the  stresses  at  the  crack 
tip  are  infinite  so  that  at  the  crack  tip 


Hence  for  a patched  crack  the  result  of  including  effects  of  the  transverse  shear  stress  in  the 
sheet  itself  is  to  reduce  the  stress  intensity  factors  by  the  factor  (1  — (tgKs/l2Gs)).  For  most 
practical  situations  this  reduction  will  only  be  of  the  order  of  a few  per  cent.  Nevertheless  this 
is  an  important  result  as  it  shows  that  the  usual  procedure  of  neglecting  the  shear  stresses  xix 
and  rzv  in  the  sheet  is  conservative  as  regards  estimates  of  the  stress  intensity  factors  and  hence 
slightly  underestimates  the  load  required  for  crack  initiation.  For  example  for  the  patch  discussed 
earlier  Kt  = 6-56  x 103  MPa/mm,  ts  = 2-286  mm,  Gs  = 27-3  x 103  MPa  and  1 — (/SA»/ 
12G»)  = 0-954,  the  reduction  in  the  stress  intensity  factor  thus  being  only  4-6  per  cent. 
Consequently  when  producing  a numerical  model  for  the  crack  patching  problem  we  can  neglect 
the  effects  of  the  transverse  shear  stresses  in  the  sheet  and  in  the  patch  with  little  loss  in  accuracy. 

The  finite  element  technique  for  analysing  the  crack-patching  problem  will  be  developed  in 
the  next  section. 

3.  A FINITE  ELEMENT  APPROACH  TO  CRACK  PATCHING 

In  the  previous  Section  we  have  seen  that  transverse  shear  stresses  txz  and  tvz  in  the  patch 
and  the  sheet  may  be  neglected;  hence  the  stiffness  matrices  for  elements  in  the  sheet  and  in  the 
patch  may  be  computed  using  standard  finite  element  routines  for  membrane-type  elements.  It  is 
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the  stiffness  matrix  of  the  adhesive  which  couples  these  two  stiffness  matrices  together.  We  will 
therefore  turn  our  attention  to  determining  the  stiffness  matrix  for  an  element  of  adhesive. 

The  strain  energy  Va  of  an  element  of  the  adhesive  is  given  by 


Va  2 Ga 


(Tzz2  + rzy2)dxdydz  = 

20a 


(tiz2  + rz/)dxdy 


(3.1) 


where  the  double  integration  is  over  the  area  of  the  element.  The  element  stiffness  matrix,  Ka'  for 
the  adhesive  is  now  obtained  by  first  prescribing  the  form  of  the  shear  stress  variation  within 
the  element  and  then  by  differentiating  the  strain  energy  with  respect  to  each  of  the  elemental 
degrees  of  freedom.  However  before  prescribing  the  nature  of  the  shear  stresses  we  must  first 
stipulate  the  geometry  of  the  element  to  be  considered.  Here  we  will  confine  our  attention  to 
triangular  elements  which  are  the  basic  elements  from  which  all  other  elements  may  be  assembled. 
Let  us  assume  that,  within  the  element,  the  shear  stresses  vary  linearly  with  x and  y,  i.e., 

Tzx  — Kfl<  + hx  + c(y) TZZi  + (a)  + bjx  + C)y)rzzl  -f  (am  + bmx  + cmy)Tzzm\llA  (3.2) 
where 


ai  — X})’m  Xm)'j 


bt=  yt  — ym 


Ci  = Xm  — X] 


(3.3) 


and  ay  am,  by  bm  and  cy  cm  are  obtained  by  a cyclic  permutation  of  the  indices  f,  j,  m.  Here 
(xi,  yi),  ( X] , yi ) and  (xm,  ym)  are  the  coordinates  of  the  corners  of  the  element,  A is  the  area  of 
the  element  and  rzzi,  tzxU  Tzzm  are  the  values  of  the  shear  stresses  at  the  nodes  /,  j and  m 
respectively  (see  Figure  3).  Hence  if  we  now  substitute  for  rzzt,  tzz,  and  Tzzm  as  given  by  equations 
(2.36)  we  obtain 


FIG.  3 ELEMENT  OF  ADHESIVE 
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similarly. 
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(/5/a  — /4 2)  LiP=t.f.m 


<t> p(x,  4’)[/t(«o  — M»)p  4-/5(i’o  — n*)p]  (3.5) 


where  <f>t(x , 4)  = (o  4-  Af.v  4-  ay)/2A  and  where  p takes  on  the  values  i,  j or  m.  Substituting 
for  txz,  and  ryz  into  equation  (3. 1 ) we  obtain  the  following  expression  for  the  strain  energy: 


I*  • 

7 = 

2Ca 

J •) 


[F]  {Zptf>p(x , y)(u0  — i/»)p}2  4-  F3  [~p  <f>p(x,  y)(v0  f«)p}2 


2F2  {-p^p(.v,4’)(i’«  — v,)p){'L Q 4>„(x, y){u0  — m»)?}]  r/Arr/y 


(3.6) 


Here  p and  q are  dummy  indices  which  take  on  the  node  values  of  either  /,  j or  m and  where, 
for  the  sake  of  convenience,  we  have  written 

Fi  = (/42  4-/e2)/(/5/a  -/t2)2 


Fz  = /4(/5  4 -/e)/(/5/e  -/42)2 


(3.7) 


Fa  = (/42  4-/52)/(/5/a  -/42)2 

Equation  (3.6)  directly  relates  the  strain  energy  of  the  adhesive  to  the  displacements  in  the  sheet 
and  the  patch  and  clearly  illustrates  the  importance  of  the  previous  Section  where  the  exact 
nature  of  the  shear  stress  in  the  adhesive  was  determined. 

The  element  stiffness  matrix  is  now  obtained  by  differentiating  Va  with  respect  to  each  of 
the  elemental  degrees  of  freedom,  i.e. 

(38) 


where  7 is  the  vector 

7 ’l'  — [w*|,  V si,  Uoi , V0iy  UhJ , Vsjy  UoJ,  Vojy  Hjinit  1'gffl,  Worn,  Pom]  (3.9) 

This  technique  is  explained  in  detail  in  references  [1,2],  and  the  full  form  of  the  12x12  stiffness 
matrix  is  given  in  Appendix  I. 

Having  thus  obtained  the  stiffness  matrix  for  the  adhesive,  we  now  couple  the  stiffness  matrix 
of  the  adhesive  to  the  stiffness  matrix  of  the  composite  patch  and  thus  produce  a stiffness  matrix 
for  the  bonded  patch-adhesive  pair.  This  coupling  is  a very  simple  procedure  since  by  specifying 
the  degrees  of  freedom  of  the  adhesive,  i.e.  7,  we  immediately  specify  the  position,  in  the  matrix, 
of  the  terms  due  to  the  composite  For  example  in  equation  (3.9)  we  see  that  the  u and  v displace- 
ments in  the  patch  occupy  positions  3,4.7, 8, 1 1,12.  Hence  the  first  row  in  the  stiffness  matrix  of 
the  patch  is  located  in  the  third  row  of  the  stiffness  matrix  for  the  bonded  element,  the  first  term 
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in  the  row  occurring  in  position  (3,3)  the  second  term  in  position  (3,4)  the  third  term  in  (3,7)  the 
fourth  term  in  (3,8)  the  fifth  in  (3,1 1)  and  the  last  term  in  (3,12).  The  second  row  is  relocated  in 
the  fourth  row,  the  third  row  is  relocated  in  the  seventh  row  etc.  An  algorithm  for  doing  this  is 
given  in  Appendix  II. 

This  bonded  pair  is  now  used  in  conjunction  with  standard  finite  element  routines  except  at 
the  crack  tip  where  it  is  coupled  to  the  special  crack  tip  element  developed  by  the  author’s  in 
references  [1,  2]. 

So  far  we  have  assumed  a complete  bond  between  the  adhesive  and  the  sheet.  However 
sometimes  one  is  faced  with  the  situation  where  debonding  of  the  patch  from  the  sheet  has 
occurred  over  some  area.  The  present  analysis  can  cover  this  case  by  allowing  all  or  any  of  t te 
nodal  displacements  at  the  lower  interface  of  the  adhesive  to  differ  from  the  displacement  in  the 
sheet.  In  the  case  of  a small  debond  on  either  side  of  the  crack  this  requires  correspondingly 
small  elements,  and  hence  a more  complex  mesh.  To  overcome  this  requirement  for  a fine  mesh 
it  is  possible  to  modify  the  bonded  element  developed  above  with  a nodal  displacement  vector  £ , 
by  forcing  the  displacement  of  the  adhesive  at  the  crack,  to  take  the  mean  value  of  the  displace- 
ment of  the  sheet  at  either  side  of  the  crack.  In  this  case  the  nodal  displacement  vector  is 

“R*  = + , Vgt+,  Uoi , V oii  Us],  Vs],  UoJ , Vo],  Usm,  Vsm,  Uom,  Vom , Ust  , Vsi  ] (3.10) 

where  the  ith  node  is  assumed  to  be  at  the  crack  and  where  Ustu  and  i/„<  etc.  are  the  values  of 

the  displacements  on  either  side  of  the  crack.  The  stiffness  matrix  Ra  for  this  case  is  computed 

using  the  algorithm 

Ka(m,  n)  — Ka(m,  n ) for  m,  n = 1 to  1 2 
Ka(m,  n)  = \Ka(m,  n ) for  m < 2,  n < 12 
Ra{m,  n)  — \Ra(m,  n)  for  m ^ 12,  n < 2 (3.1 1) 

£„(/n,  12  -f  n)  = Ra(.m,  n)  for  n < 2,  m ^ 12 

Ra(m  + 12,  n)  = Rtt(m,  n ) for  m < 2,  n < 14 

where  Ka(m,  n)  is  the  m,  nth  term  of  the  stiffness  matrix  determined  above  assuming  complete 
adhesion. 

It  is  interesting  to  note  that  in  the  case  of  a centre  cracked  panel  under  uniform  tension  this 
results  in  the  adhesive  at  the  crack  h iving  zero  displacement,  whereas  in  the  fully  bonded  situation 
the  adhesive  will  move  differently  on  either  side  of  the  crack.  The  latter  results  in  a displacement 
discontinuity  in  the  adhesive  across  the  crack  which,  in  turn,  may  result  in  the  creation  of  a crack 
in  the  adhesive  itself  or  may  cause  a localized  debonding  effect  along  the  crack. 

Having  thus  allowed  for  all  possible  effects  we  will  now  illustrate  the  use  of  the  analytical 
tools  developed  above  by  considering  a tension  panel  with  a central  crack  which  is  reinforced  by 
a Boron-epoxy  patch  and  where  we  will  concentrate  on  the  reduction  in  the  stress  intensity  factor 
caused  by  the  patch;  and  on  the  stresses  in  the  patch  and  the  adhesive. 


4.  ILLUSTRATION  OF  THE  METHOD 


As  mentioned  in  the  last  paragraph  above,  let  us  consider  a thin  rectangular  sheet  of 
dimensions  508  mm  x 635  mm  x 2-3  mm  subjected  to  a uniform  tensile  stress  of  689  KPa. 
The  sheet  contains  a central  crack  38  - 1 mm  long  which  is  patched  on  both  sides  and  at  both  ends 
by  a Boron-epoxy  laminate  placed  directly  over  the  crack  tip.  The  laminate  is  uniaxial,  the  fibers 
lying  perpendicular  to  the  length  of  the  crack,  the  lateral  dimensions  of  the  patch  are  100-8  mm  x 
17-7  mm,  and  various  thicknesses  of  patch  are  considered.  The  patch  is  bonded  to  the  sheet  with 
an  adhesive  of  thickness  0- 1016.  The  Young's  modulus  of  the  sheet  is  71  -02  x 103  MPa  while 
the  Poisson’s  ratio  is  0-32.  The  moduli  of  the  laminate  are  G'23  = 4-94  x I03  MPa,  G is  = 
7-24  x 103  MPa,  £j  = 208-1  X 103  MPa,  Ei/Ez  = 8-18,  vi2  = 0-1677  and  the  shear  modulus 
of  the  adhesive  is  Ga  = 9-65  x I02  MPa. 

Because  of  the  symmetrical  nature  of  the  problem  only  one  half  of  the  structure  was  analysed 
and  the  sheet  was  idealised  as  in  figure  4 where  in  order  to  model  accurately  the  stresses 
in  the  adhesive  it  was  necessary  to  have  a fairly  fine  mesh.  Indeed,  by  modelling  one  half  of  the  plate 


rather  than  only  one  quarter  we  are  able  to  calculate  both  the  mode  I and  mode  II  stress  intensity 
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factors,  Ki  and  Aj,  for  the  structure,  as  described  in  references  [1,2]  so  that  any  error  in  the  data 
can  be  detected  by  evaluating  the  ratio  K-i/K  1.  If  the  data  are  correct,  then  this  ratio  should  be 
zero.  However,  because  of  the  large  amount  of  data  and  because  of  round-off  effects,  which 
accumulate  during  computation  it  was  often  found  that  the  ratio  was  non  zero.  The  maximum 
tolerated  value  of  this  ratio  was  0 01,  i.e.  K-i  was  not  allowed  to  exceed  1%  of  K\.  Whenever  A* 
did  exceed  this  value  it  transpired  there  were  data  errors. 

Table  I shows  the  effect  that  increasing  the  patch  thickness  has  upon  the  computed  value  of 
the  stress  intensity  factor  Ai,  when  no  debond  is  present. 

TABLE  1 


Patch  Thickness  (mm) 

Ai  (KPa.m1) 

0 

166-2 

0-127 

77-8 

0-254 

65-0 

0-381 

57-9 

0-508 

53- 1 

0-762 

46-7 

It  is  interesting  to  note  that,  when  the  patch  is  not  present,  the  method  yields  the  value  Ai  = 
166-2  kPa.m1  which  differs  by  only  1%  from  the  exact  analytical  result6  168-5  kPa.m*.  From 
Table  1,  we  see  that  patching  is  a very  efficient  way  of  lowering  the  stress  intensity  factor, 
and  that  the  addition  of  extra  layers  of  the  Boron-epoxy  laminate  continually  decreases  the 
value  of  the  stress  intensity  factor.  However,  increasing’  the  thickness  of  the  laminate  adversely 
effects  the  Tty  shear  stress  in  the  adhesive  as  may  be  seen  in  Table  2 where  the  maximum  shear 
stress  developed  in  the  adhesive  at  the  ends  of  the  patch,  is  shown  for  various  patch  thicknesses. 


TABLE  2 


Patch  Thickness  (mm) 

Tzy  (kPa) 

012 

106 

0-254 

181 

0-371 

236 

0-508 

280 

0-762 

345 

From  this  table  it  is  apparent  that  increasing  the  thickness  of  the  laminate  from  0-5  mm  to 
0-76  mm  (i.e.  from  four  to  six  layers)  results  in  an  increase  in  the  shear  stress  in  the  adhesive  of 
23%  while  only  reducing  A,  by  12%. 

Hence,  when  increasing  the  number  of  layers  of  the  laminate  one  must  be  watchful  of  the 
detrimental  effect  it  has  in  raising  the  shear  stress  in  the  adhesive.  One  must  therefore  be 
careful  when  using  composite  patches  not  to  exceed  the  “optimum"  number  of  layers. 

Figures  5,  6,  7 and  8 show  the  complete  variation  of  the  shear  stress  throughout  the  adhesive 
for  various  thicknesses  of  the  laminate.  From  these  diagrams  it  is  apparent  that  the  shear  stresses 
in  the  adhesive  on  the  line  of  the  crack  at  the  right  hand  edge  of  the  patch  (i.e.  outboard  of  the 
crack  tip)  are  much  smaller  than  those  at  the  left  hand  edge  of  the  patch  (i.e.  inboard  of  the 
crack  tip).  Also,  the  stresses  at  the  right  hand  edge  vary  along  the  length  of  the  patch  in  the  same 
general  fashion  as  in  a lap  joint. 


6.  G.  C.  Sih  and  Mathematical  Fundamentals,  Fracture,  Vol.  2,  Academic  Press,  N.Y., 

H.  Liebowitz  1969. 


12 


) 

l 


Left  hand  side  of  patch 


FIG.  5 VARIATION  OF  ADHESIVE  SHEAR  STRESS:  THICKNESS  OF  PANEL  = 0.012/  mm 


Left  hand  side  of  patch 


G.  6 VARIATION  OF  ADHESIVE  SHEAR  STRESS:  THICKNESS  OF  PANEL  = 0.2540  mm 


Middle  Patch 


VARIATION  OF  ADHESIVE  SHEAR  STRESS:  THICKNESS  OF  PANEL  = 0.5080  mm 


The  stresses  in  the  patch  are  predominantly  ay  stresses,  i.e.  ax  « TIy  as  0,  and  are  fairly 
constant  across  the  width  of  the  patch,  the  only  exception  being  at  the  crack,  i.e.  on  the  x axis. 
These  stresses  are  shown  in  Figure  9 where  we  see  that  for  each  thickness  of  patch  the  maximum 
stresses  occur  in  the  region  of  the  crack  and  decay  monotonically  towards  the  edges  of  the  patch. 
The  maximum  ay  stress  is  shown  in  Table  3 for  various  thicknesses  of  the  patch. 


TABLE  3 


Patch  Thickness  (mm) 

(kPa) 

0127 

10490 

0-254 

6780 

0-381 

5128 

0-508 

4159 

0-762 

3040 

Knowing  the  stress  distribution  in  the  patch  it  is  possible  to  calculate  the  load  carried  by  the 
patch.  In  this  case  it  is  of  interest  to  note  that  increasing  the  thickness  of  the  patch  from  0-5  mm 
to  0-76  mm  increases  the  load  in  the  patch  by  only  9 -6%. 

Let  us  now  examine  the  effect  of  a debond  between  the  plate  and  the  adhesive;  the  debond 
is  assumed  to  occur  only  in  the  vicinity  of  the  crack  and  has  a width  of  5 08  mm.  For  each  patch 
configuration  mentioned  above  the  effect  of  the  debond  on  the  stress  intensity  factor  and  on  the 
stresses  in  the  patch  was  less  than  0-1%.  The  major  effect  of  the  debond  was  to  lower  the  rty 
stresses  in  the  adhesive  at  the  crack.  This  effect  can  be  seen  in  Figures  5,  6,  7 and  8.  However,  as 
can  also  be  seen  in  these  Figures,  the  shear  stress  at  the  edge  of  the  debonded  region  is  virtually 
unaffected  by  the  existence  of  a debond. 

These  comments  apply  only  to  a relatively  small  debond  along  the  crack;  however,  there  is 
no  physical  reason  to  suppose  that  a debond  does  in  fact  exist,  indeed  as  mentioned  in  reference 
[3]  the  presence  of  a debond  around  the  crack  in  an  actual  experimental  specimen  has  not  yet 
been  reported.  Consequently  if  a debond  does  exist  its  width  must  be  small.  The  theoretical 
existence  of  a debond  is  totally  dependent  upon  the  theoretically  zero  width  of  a crack  and  the 
resultant  stress  discontinuity,  in  the  adhesive,  across  the  crack.  However,  in  practice  the  crack 
has  a small  but  finite  width  so  that  this  discontinuity  will  not  occur.  The  shear  stress  in  the 
adhesive  will  still  be  high  in  the  vicinity  of  the  crack  and  in  the  normal  range  of  the  working 
stresses  the  adhesive  may  go  plastic.  Indeed  since  Figures  5,  6,  7 and  8 show  that  the  shear  stress 
in  the  adhesive  at  the  crack  is  of  the  same  order  of  magnitude  as  the  shear  stresses  at  the  edge  of 
the  patch  the  authors  believe  that,  as  discussed  by  Hart  Smith7,  in  the  normal  range  of  the 
working  stresses  the  adhesive  close  to  the  ends  of  the  patch  and  the  adhesive  directiy  over  the 
crack  may  be  behaving  plastically,  and  that  debonding  may  not  occur. 

As  a result  of  this  analysis,  the  authors  believe  that  failure  of  a patched  crack  will  primarily 
be  determined  by  the  same  failure  modes  as  reported  in  reference  [7],  These  are  as  follows. 
First,  as  a result  of  the  adhesive  going  plastic  the  rate  of  increase  of  the  stress  intensity  factor 
with  load  will  be  dramatically  increased  so  that  the  applied  load  may  cause  A,  to  exceed  Air,  the 
fracture  toughness  of  the  plate.  Second,  the  ultimate  shear  strain  of  the  adhesive  may  be  exceeded 
resulting  in  bond  failure  (this  was  the  usual  experimental  failure  mode  reported  by  Hart  Smith 
[7]).  Finally,  the  peel  stresses  induced  at  the  ends  of  the  patches  may  cause  failure,  or  the  com- 
posite patch  may  itself  fail. 


7.  L.  J.  Hart  Smith  Analysis  and  design  of  advanced  composite  bonded  joints,  NASA 

CR-22I8,  August  1974. 
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/Applied  stress 


FIG.  9 VARIATION  OF  PATCH  STRESS  WITH  DISTANCE  FROM  CRACK  (Left  hand  side  of  patch) 


APPENDIX  I 


Determination  of  the  Stiffness  Matrix  for  the  Adhesive  Element 


From  equation  (3.6) 


F«  — 2q  [Fi(i-P <f>p(u0  ^p(t’o  — t>*)p)2 

~F  2F2(lp<^p(Uo  Ux)p){^q  <f)q(Vo  Vs)q)]dxdy  (i) 

The  first  row  in  the  stiffness  matrix  is  obtained  by  differentiating  with  respect  to  uKi  (i.e.  usp 
with  p i ) which  represents  the  value  of  ux  at  the  /th  node.  Differentiating  we  obtain 


r f*  I* 

^ a ta 

hi  Ga  , 


t F i <f>p(ux  — Mo)p)  4"  F ‘Z^ii^q  <f>q(V8  — Vo))  dxdy  — 2-n  K\n^n  (ii) 


where  A„  is  the  nth  term  in  the  vector  A as  given  by  equation  (3.9)  and  K\n  is  the  nth  term  in  the 
first  row  of  the  stiffness  matrix  of  the  adhesive.  Following  the  general  technique  the  stiffness 
matrix  K is  found  to  be  as  shown  in  Table  I A. 

Here,  for  the  sake  of  convenience,  we  have  written 

< /',  / > = if  <f>i<l>)dxdy  (iii) 


</,  m> 


= <j>i<l>mdxdy 


where  the  functions  <fri,  <j>)  and  <f>m  are  as  given  by  equation  (3.5).  Following  this  general  technique 
the  stiffness  matrix  Kea  is  found  to  be  as  shown  below: 


APPENDIX  II 


Stiffness  Matrix  for  Combined  Patch-Adhesive  Element 

Define  a vector  8(/'),  where  i takes  the  values  I to  12,  such  that  8(3)  I,  8(4)  2,  8(7)  = 3, 

8(8)  4,  8(  1 1 ) 5,  8(  12)  = 6 and  8 (i)  is  zero  for  all  other  values  of  /.  Also  define  a vector  R(i) 

such  that  R(i)  is  zero  whenever  3(/)  is  zero  and  /?(/)  = 1 whenever  8 (/)  is  not  zero.  Then  if  the 
stiffness  matrix  of  the  adhesive  is  denoted  by  Ka,  and  the  stiffness  matrix  of  the  patch  is  denoted 
by  Kp,  the  stiffness  matrix  of  the  patch-adhesive  pair  K is  given  by 

K(i,i)  = KaOJ)  + R(i)R(j)KP(S(i ),  S(j))  (i) 

where  / and  j take  on  the  values  I to  12. 
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